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Abstract — In this paper we report our systematic study of 
a promising absolute calibration technique of analog photo- 
detectors, based on the properties of parametric down conversion. 
Our formal results and a preliminary uncertainty analysis show 
that the proposed method can be effectively developed with 
interesting applications to metrology. 

Index Terms — Detectors, calibration, optical parametric am- 
plifier, optical signal detection, optics, quantum theory 

I. Introduction 

In optical radiometry primary standards are based on ab- 
solute sources or detectors [1]. Synchrotron and blackbody 
radiators are absolute sources. The relative uncertainty of 
both these sources is about 1 part in 10 3 . Among the ab- 
solute detectors, there exist the following two types: electrical 
substitution radiometers (ESR) based on thermal effects, and 
semiconductor photon detectors. Uncertainties of a few parts in 
10 4 appear to be the limit of these techniques. However, such 
a high accuracy of traditional calibration methods is reached 
only for an appropriate choice of wavelength and special types 
of detectors. 

Recently the use of photons produced by means of spon- 
taneous parametric down-conversion (SPDC), where photons 
are emitted in pairs strongly correlated in direction, wave- 
length and polarization, has been exploited for the absolute 
calibration of detectors in photon-counting mode [2]-[7]. 
This absolute technique is becoming attractive for national 
metrology institutes to realize absolute radiometric standards. 
In fact it relies simply on the counting of events, involves a 
remarkably small number of measured quantities, and does not 
require any reference standards. 

Because of the success of the SPDC scheme for calibrating 
single-photon detectors, it is important to analyze the pos- 
sibility to extend this technique to higher photon fluxes for 
calibrating analog detectors. A seminal attempt in this sense 
was made in [8] following the theoretical proposal of [3]. 
Nevertheless, these results were limited to the case of very 
low photon flux (as we will show later). 

In this paper we report our systematic analysis of the 
measurement method for increasing photon flux produced 
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by means of SPDC, showing how to estimate the quantum 
efficiency in a real analog regime. Furthermore, we point 
out the intrinsic limitation of SPDC for calibration when 
one moves towards fluxes requiring large parametric gain. 
Thus, the stimulated parametric down conversion (PDC) is 
considered as an alternative bright source of correlated beams 
and we demonstrate that, under some opportune conditions, 
the quantum efficiency can be estimated in this new regime. 

II. The absolute calibration technique 

The scheme for calibrating photon detectors by using para- 
metric down conversion is schematically depicted in Fig. 1. 
It is based on the specific properties of this process, where 
a photon of the pump beam (usually a laser beam) "decays" 
inside a non-linear crystal into two lower-frequency photons, 
1 and 2 , such that energy and momentum are conserved 
(Wp Um! , = cJi + u)2, k pump = k\ + &2). Moreover, the two 
photons are emitted within a coherence time r co h of tens of 
femtoseconds from each other. The process can be spontaneous 
(SPDC) when no modes of radiation except the pump modes 
are injected through the input face of the crystal. If a seed 
mode &2 is injected, its presence stimulates the process and 
many more photons of the pump are converted. 

In essence, the calibration procedure consists of placing 
a couple of detectors Dl and D2 down-stream from the 
nonlinear crystal, along the directions of propagation of cor- 
related photon pairs for a selected pair of frequencies. Since 
the incident photon fluxes Fi(t) and i 7 ^) are correlated 
within 1CP 13 s, the fluctuations of the registered currents 
ii(t) and (i) are suppose to be strictly correlated. The 
non ideal quantum efficiency of the detectors makes some 
photons missed sometimes by Dl sometimes by D2, spoiling 
the correlation. The techniques for estimating the quantum 
efficiency, both in counting and in analog regime, consist in 
measuring this effect. 

In the following the photodetection process in the analog 
regime will be modelled as a random pulse train [9] 

i(t) = ^2q n f(t-t n ), 

n 

i.e. a very large number of discrete events at random times 
of occurrence t n . The pulse shape /(<) is determined by 
the transit time of charge carriers. We assume that f(t) is 
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Fig. 1. Scheme for calibration of light detectors. The processing of the 
detectors output indicated here is relative to the analog regime. 



a fixed function with the characteristic width r p and a unit 
area (f(t) has dimension s _1 ). t p ~ 10 ns represents a 
typical value in analog detection. The pulse amplitude q n is 
a random variable representing the electron charge produced 
in the n-th event, in order to account for a possible current 
gain by avalanche multiplication. The statistical nature of the 
multiplication process gives an additional contribution to the 
current fluctuations [10]. In an ideal instantaneous photocell 
response, without avalanche gain, all values q n are equal to 
the charge e of a single electron and f(t) ~ S(t). In the case 
of ideal quantum efficiency, since the probability density of 
photo-detection at time t at detector Dk (k — 1, 2) is related 
to the quantum mean value (F k (t)) of the photon flux (photons 
per second), we calculate the average current output of Dk as 
(see for example [11]) 

(i k ) = J2(lknf(t - t kn )) = I dt k (q k )f(t - t k ){F k (t k )) 
n J 

(1) 

where the factor (q k ) is the average charge produced in a 
detection event. We have assumed the response function for 
the two detectors to be the same, /i(t) = fe(t) = f(t). 

Now we introduce the quantum efficiency r) k of detector 
Dk, defined as the number of pulses generated per incident 
photon. In [12] a real detector is modelled with an ideal 
one (ij = 1) preceded by a beam splitter with transmission 
coefficient equal to the quantum efficiency of the real detector 
[13]. Following the results reported there we can perform 
in Eq. (HJ the substitution F k (t) — ► rj k F k (t). Thus, being 
(F k (t)) time independent, we obtain: 



respectively for k = j where (F k (t k )F k (t' k )) is the auto- 
correlation function of the photon flux at detector k, and for 
k 7^ j where (F k (t k )Fj(tj)) is the cross-correlation between 
the fluxes incident at the two different detectors. 

A. Analog calibration using SPDC 

The functions (F k (t k )F k (t' k )) and (Fk(tk)Fj(tj)) are eval- 
uated in [12] for SPDC. As mentioned above, the resolving 
time of a real analog detector is finite, and in general can be 
considered much larger than the SPDC coherence time. Thus 
any fluctuations in the intensity of light are integrated over t p 
during the detection process. In the limit of r p 3> r co ^, for 
k = 1 and j = 2 we have 



(h(t)h(t + r)) = (t,) 2 + m (qf)Hr) ■ [(Fx) ^ 

(ii(t)i 2 (t + t)> = (ii)(ia) + ^2(91) (^(r) 



(4) 



(5) 



(ik) = Vk(qk)(F k ) . 



where we have introduced the convolution of the response 
function of detectors J-{t) = J dtf(t)f(t + r). The term 9 
depends on the second power of the parametric gain V, i.e. the 
mean number of photons per mode of the radiation. For our 
purposes we can estimate approximatively 3 ~ V(F) [12]. 
The Eqs (0]i and <(5j are the fundamental tools for studying the 
problem of absolute calibration of analog detectors and thus 
we are going to discuss them in detail. Despite they seem to 
be quite symmetric, we observe an important difference. The 
term proportional to (Fi ) in the autocorrelation function is the 
shot noise contribution and for this reason the quantum effi- 
ciency 7/1 enters linearly, while in the current cross-correlation 
function the corresponding term is due to the high quantum 
correlation between the signal and idler beams of PDC and 
the quantum efficiency appears quadratically. Its presence is 
the key for absolute calibration. The term 3 both for auto- 
and cross-correlation functions, can be neglected as long as 
V <C 1, i.e., the mean number of photons per mode is much 
smaller than one. However, if the duration of the photocurrent 
pulse is much longer than the coherence time, this assumption 
does not prevent photodetection being in a strongly analog 
regime, because a lot of photons can be absorbed during the 
pulse duration generating overlapping of pulses [see Eq (HJ]. 
The term proportional to (i) 2 oc (F) 2 is due to the presence 
of more than one photon within a time interval r p and for that 
reason is more delicate. It can be neglected in Eqs. (0]l and (0 
only if (F) <C F(t). Since the pulse f(t) has a height around 



(2) 1/t p 



[^(r)] = .F(O) ~ 1/t p . So the condition becomes 



The auto-correlation and the cross-correlation functions for the 
currents can be expressed as 



(F)t p <C 1, i.e., the number of incident photons during the 
resolving time of the detector should be much less than one, 
i.e. one should work in a non-overlapping regime. 
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Therefore, it is convenient to distinguish among three dif- 
ferent regimes: 

(I) (F)t p -C 1 (i-e. photocurrent pulses do not overlap on 
the average). For a detector with a time constant r p = 10 ns, 
the corresponding photon flux must be below 10 s photons/s, 
which, for a wavelength of 500 nm, means a power of tens of 
pW. From and (0, neglecting 3 and the terms in (i) 2 one 
obtains a formula for the estimation of quantum efficiency like 
the one reported in [3]. However, this regime does not appear 
interesting because of the necessity to work only at very low 
intensity, where no overlapping between pulses happens. In 
principle, one could distinguish between different pulses of the 
current and work in the counting mode provided the amplitude 
q n of each pulse is large enough to be detectable. 

(II) (F)t p > 1 but still V <€. 1 (i.e., photocurrent pulses 
overlap but the parametric gain and photon flux are still quite 
low). By considering the same parameters as used in case I, 
coherence time r co /j of the order of 100 fs, and the requirement 
that V < 0.001, this means a photon flux up to 10 10 photons/s, 
i.e. a power of few nW. Here, only the term 3 can be neglected 
in Eqs. and (0. Defining the current functions as 8ik{t) = 
ik{t) — (ik), the analog quantum efficiency can be estimated 
as: 

(qf) s {Sh(t)6i 2 (t + T)) 
mm) = 7— (qi ) /x . Mx . ,, , — rr. (6) 
WIT (foi(f)<hi(t + T)) 

Unfortunately it is not satisfying from the metrology point 
of view because of the presence of an unknown parameter 
M = (q 2 )/(qi} 2 related to the statistics of charge fluctuations 
and that has to be estimated independently in some other 
way. We suggest to avoid the problem by integrating Eq. (0 
over time r. It correspond to evaluate the power spectrum 
of the fluctuations at frequencies around zero, namely much 
smaller than 1/t p . We would like to stress that in this case, 
the assumption /i(t) = f2(t) = f(t) is not necessary. Since 
J (1tJ-(t) = 1, integrating Eq. ® in r and dividing by the 
current [see Eq (fJJ], we obtain 

/ dT(6ii(t)5i 2 {t + T)) 



(7) 



Here r\ 2 (92} is the electron charge produced per single incident 
photon, or, according to formula (0, the ratio between the 
current and the photon flux. 

Eq. (0 represents one of the main results of the paper, since 
it shows that the absolute calibration of analog detectors by 
using SPDC is indeed possible. 

(Ill) V > 1 (fluxes larger than 10 10 photons/s). In this 
regime each term of (0 and (0 is important and no simple way 
can be found for the absolute calibration with SPDC. It can be 
shown that in this case one should be able to collect exactly 
the same number of correlated modes by Dl and D2 [12]. 
Since SPDC takes place with a very large spectral and spatial 
bandwidth, it would require accurate and well balanced spatial 
and frequency selection. This could originate systematic errors 
which are difficult to be evaluated. 

B. Analog calibration using stimulated PDC 

We overcome the problem of calibration for photon fluxes 
larger than 10 10 , highlighted in point (III), exploring the 



possibility of using a parametric amplifier configuration, where 
a seed coherent beam with power (f> is injected along direction 
2 stimulating the bright emission of the two correlated beams. 
In this case the photon fluxes can be increased by varying 
the power <f> of the coherent seed beam, without increasing 
V (we can lead back to the condition V <C 1). Introducing 
the proper expression for the fluxes in the stimulated emission 
in (0, under the condition that the spontaneous emission is 
negligible we have 



(«i)s = m(qi)v<f> 

= m(q2){i + v)4> 



(8) 



where the subscript s reminds that here the currents are cal- 
culated for the stimulated PDC. Using Eq. the correlation 
functions of the current fluctuations can be expressed in a 
simple form, under the assumption of V <C 1: 

(ShiQSitit + r)). = m (ql)T(T)V<p 

(5i 2 (t)5i 2 (t + T)) s = m (ql)F{ T ){l + V)d> (9) 

(Sh(t)Si 2 (t + T)) s = 2 mm { qi ){q 2 )T(T)V<P 

The quantum efficiency can be evaluated as 



, N 1 (q 2 i) 1 , {Sh(t)Si 2 (t + T)) B 

2{qi) z (5ii(t)5ii{t + T)) 8 



(10) 



Alternatively, integrating in r the expression for the cross- 
correlation we have 



. 1/ dr<fli(t)fi 2 (t + T)) B 

m{Q2) = ■x 7-t 

2 n s 



(11) 



in which the avalanche gain factor disappears. There are 
two main advantages of using the stimulated PDC. First of 
all we are not up-limited in the photon fluxes that we can 
use. Considering detector without internal gain where the 
Johnson noise in the amplification process is predominant, it is 
important to obtain a good signal to noise ratio by increasing 
the signal. We stress that in our case the signal is of the order 
of the shot noise! A second advantage lies in the very narrow 
bandwidth, both in space and frequency of the stimulated 
emission, that is the one of the coherent seed beam. This makes 
easier to collect the correlated beams resulting in two bright 
spots in the detection plane. 

III. Preliminary uncertainty analysis 

The calculation of the expected uncertainties in the evalua- 
tion of the correlations functions depends on the measurement 
technique. Here we consider a measurement process in which 
the two current outputs of Dl and D2 are combined by an 
analog multiplier and the result is integrated for a time T much 
larger than the response time t p . The general calculation of 
the variance of this quantity would require the knowledge of 
the fourth-order quantum correlation functions of the photon 
fluxes which are rather complex. Therefore we used the 
approximation of small fluctuations, i.e. Si(t) <C (i), which 
perfectly fits the regime of strong overlapping discussed in 
the previous section. In this case one obtains, both for the 
auto- and cross-correlation, a maximal uncertainty 
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A^^^^T- 1 / 2 (12) 

Considering equation (|6]l for a detector without avalanche 
gain, following the uncertainty propagation rule, we obtain the 
relative uncertainty on the quantum efficiency as 




where iV r is the number of photons detected during the 
detector response time. The uncertainty scales with the square 
root of the measurement time T. For example, if T = 1 s, it is 
lower than one part in 1CP 3 . Since, as we are going to discuss, 
this term is not the dominant one this estimate suffices for our 
purposes. 

The noise contributions from the detector (dark current, 
noise of the transimpedance amplifier) and background light 
are supposed to be statistically independent of each other 
and of the photocurrents produced by SPDC light in the two 
detectors. Thus, they do not give any contribution to the cross 
correlation function of the two current fluctuations in the 
numerator of ©. In the denominator of the same equation, 
containing the autocorrelation of the current fluctuations, their 
effect should be measured (for example, rotating by 90° the 
polarization of the laser pump, to eliminate the SPDC signal), 
and subtracted. The statistical uncertainty of this measurement 
can be kept so small that it does not increase significantly the 
total relative uncertainty given above. 

The main systematic contribution is due to the non linear 
crystal optical losses. The uncertainty in its measurement 
could be reduced to few parts in 10 3 [14], [15]. Thus, it is 
the dominant contribution to the uncertainty budget, largely 
exceeding the statistical one. 
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IV. Conclusion 

Motivated by the necessity of a general absolute calibration 
scheme for detectors for various applications, we have per- 
formed a systematic study on the possibility of applying PDC 
calibration methods to the analog regime. Our results show 
that the measurement of correlations between detector output 
currents can indeed be used to extend the absolute calibration 
method to the analog regime. In particular, it is shown that 
integration of the photocurrent correlation functions in time 
allows one to avoid the measurement of the photocurrent pulse 
shape and to eliminate the necessity to know the statistics 
of the internal gain of the detector. Also, our analysis shows 
that is possible to go beyond the regime of non-overlapping 
photocurrent pulses, which was used in earlier works, and 
to move on to higher intensities. Furthermore, the stimulated 
operating mode of PDC has been considered with promising 
perspectives. 

A preliminary theoretical uncertainty analysis gives a rela- 
tive value around 10~ 3 , mainly due to the measurement un- 
certainty of the optical loss in the crystal, which is interesting 
from a metrology view-point. 



